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ABSTRACT 

We report on the discovery of several compact regions of mid-infrared emission in 
the starforming circum nuclear disk of the starburst/Seyfert 2 galaxy NGC 7582. The 
compact sources do not have counterparts in the optical and near-infrared, suggesting 
that they are deeply embedded in dust. We use the [Nell]12.8u.m line emission to 
estimate the emission measure of the ionized gas, which in turn is used to assess the 
number of ionizing photons. Two of the brighter sources are found to have ionizing 
fluxes of ~ 2.5 x 10 52 , whereas the fainter ones have ~ 1 x 10 52 photons s _1 . Comparing 
with a one Myr old starburst, we derive stellar masses in the range 3-5 xlO 5 M and 
find that the number of O-stars in each compact source is typically 0.6-1. 6xl0 3 . We 
conclude that the compact mid-infrared sources are likely to be young, embedded star 
clusters, of which only a few are known so far. Our observation highlights the need for 
high resolution mid-infrared imaging to discover and study embedded star clusters in 
the proximity of active galactic nuclei. 

Key words: Galaxies: individual: NGC 7582; Galaxies: star clusters; Infrared: galax- 
ies; Galaxies: Seyfert 
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• • . 1 INTRODUCTION 



Starburst activity in nearby galaxies has, with the supe- 
rior resolution of the Hubble Space Telescope (HST), been 
resolved into disc rete massive young clu sters since the pi- 
oneering work by iHoltzman et all il992T) . In the literature, 
these clusters are usually referred to as super star clusters 
(SSC) or young massive clusters (YMCs). Their masses are 
typically ^ 10 5 M , radii < 5 pc and ages < 100 Myr. Most 
known you ng clusters in the Milky Way have masses below 
- 10 4 M jFieer et alJll999h . but there also exist examples 
of YMCs in the Milky Way, such as the Westerlund 1 clus- 
ter which has a mass of ~ 10 s M llClark et al]l2005T) . The 
most nearby analog to a YMC is the 30 Dor cluster in the 
Large Magellanic Cloud. 

YMCs are blue compact sources, whose light is dom- 
inated by massive stars. As in the case of massive star 
association formation, the progenitors of YMCs are ther- 
mal radio sources deeply embedded in dust, e.g. the ul- 
tra dense Hn regions, first discovered in Henize 2-10 by 
iKobulnickv fc Johnson! (|l999). These objects are sources of 
mid-infrared (MIR) continuum and bright nebular line emis- 
sion. Because of high dust extinction, they are usually very 
faint in the visible, and even in the near-infrared (NIR). 



* Based on observations obtained with the ESO Very Large Tele- 
scope 
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The first evidence for the existence of such MIR sources was 
unveiled by the Infrared Space Observatory (ISO) in the 
Antennae Galaxies NGC4038/9 jMirabel et al.lll998ft . Most 



of the MIR flux from this galaxy complex comes from a 
small, optic ally insignifica n t, reg ion. Subsequent NIR spec- 
troscopy bv lGilbert et alJ i200(J) suggests that the related 
cluster formed about 4 Myr ago, has a mass of «16xl0 6 
Mq, contains massive stars with T e ff up to 3.9 x 10 4 K, and 
is deeply embedded in dust, with A v ~ 10. 



Embedded YMCs have so far been observed in only 
a few galaxies: the A ntennae and Henize 2-10, as a lready 
mentioned, NGC 5253 JGoriian. Turner fc Beckl200ll). SBS 
0335- 052 JPlante fc Sauvaed l2002p IIZw40 iBeck et alJ 
12002ft and NGC1808 and NGC1365 jGalliano et al.ll2005ft . 



Here we report on the discovery of bright, compact 
MIR sources in the composite starburst/Seyfert2 galaxy 
NGC 7582, using a VLT/VISIR image centered on the 
[Nen]12.8u.m emission line. We present the data in terms 
of suggestive evidence that the MIR sources are embedded 
YMCs, and use the [Neil] luminosity to estimate physical 
properties, such as the number of ionizing photons and the 
number of massive stars within each source. We also em- 
phasize that MIR imaging at high resolution is mandatory 
to study embedded star formation in the vicinity of active 
galactic nuclei (AGN). 
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2 THE DATA 

The data were obtained at the VLT-Melipal telescope with 
VISIR, a combined imager and spectrograph for observa- 
tions in the N and Q bands (corresponding to the two at- 
mospheric windows at 8-13 and 16.5-24.5u.rn, respectively). 
A high-resolution spectrum and one narrow-band image cen- 
tered on the [Nen]12.8u.m line were obtained during the pe- 
riod September 29 to October 2, 2004. The VISIR detector 
is a SiAs 256x256 DRS array, and a pixel scale of 0'.'127 
was used, yielding a field of view of 32'.' 5x32'.' 5. Total in- 
tegration time was one hour for both image and spectrum, 
and conditions on both nights were clear with humidity be- 
low 10%. The spectrum was taken in the high-resolution 
[Nen]12.8u.m Echelle mode using a 0'.'75 wide slit, giving a 
spectral resolution R « 16000, or ~22 kms" 1 . The diffrac- 
tion limit is reached at 12.8u.rn, hence both the spectrum 
and the image have the optimal spatial resolution of 0'.'38. 
More details regarding the data are presented by Wold et 
al. (in preparation). 

NGC 7582 has a redshift of z = 0.00525, hence one 
arcsec corresponds to w 108 pc in the galaxy rest frame 1 , 
and the [Nell] line is redshifted to 12.877u.rn. The [Nell] 
narrow-band filter has A c = 12.81 u.m and half band width 
0.21 u.m, thus encompassing [Nell] emission from the galaxy. 

The [Nell] narrow-band image displayed in Fig.Qshows 
the central 6"x6" around the AGN. The MIR emission is 
resolved into several components, labelled M0-M6, where 
M0 corresponds to the AGN. There are two noticeably 
bright, dense knots of [Nell] emission to the South, labeled 
Ml and M2. The central AGN is unresolved whereas the full 
width at half maximum (FWHM) of Ml and M2 are 0'.'60 
and 0'.'65, respectively. M2 appears slightly extended in the 
SE-NW direction. The compact sources are themselves em- 
bedded in a larger and more diffuse envelope roughly 300- 
400 pc across. The position of the slit is indicated in the 
image, and the high-resolution spectrum of the [Nell] line 
from Ml is shown to the left. 

The measured fluxes and flux densities of the compact 
sources are listed in Table Q The aperture diameters used 
correspond to lx and 2x the diffraction- limited resolution. 
We estimated [Ne n] line fluxes from the narrow-band image 
by calibrating the image using the Ml line flux from the 
spectrum. The accuracy of the line fluxes derived from the 
narrow-band image in this manner is approximately 10 %. 
The spectrum shows that there is no continuum detected 
from Ml. The noise level is « 0.09 Jy, so the 3<r upper limit 
on the continuum is 0.27 Jy. By calibrating the narrow-band 
image using the Ml spectrum, we assume that the other 
sources M2-M6 also have non-detectable continua. 



3 COMPARISON WITH DATA AT OTHER 
WAVELENGTHS 

NGC 7582 is known to have a circum nuclear kpc-scale 
disk with active star formation, containing significant 
amounts of dust iMorris et alJ ll98St iReean fc Mulchaevl 
Il999l: ISosa-Brito et alJl200ll) . This is however the first time 
that the disk has been viewed at high resolution in the MIR, 

1 We assume Hq = 70 kms -1 Mpc" 1 



Table 1. Measured quantities. Columns (2) and (3) show flux 
density within aperture diameters 0''38 and 0''76, respectively. 
Columns (4) and (5) show line flux, also within apertures of 0'.'38 
and Of.' 76. 



Source [Nell] flux density [Nell] flux 

(mjy) (lCT 16 Wm" 2 ) 



(1) 


(2) 


(3) 


(4) 


(5) 


Ml 


45.3±4.5 


131.2±7.7 


1.5±0.1* 


2.5±0.1* 


Ml 


45.3±4.5 


131.2±7.7 


1.50±0.15 


4.35±0.26 


M2 


43.9±4.5 


126.7±7.6 


1.46±0.15 


4.20±0.25 


M3 


29.4±3.6 


103.2±6.8 


0.98±0.12 


3.42±0.23 


M4 


29.0±3.6 


95.9±6.6 


0.96±0.12 


3.18±0.22 


M5 


23.1±3.2 


68.3±5.6 


0.77±0.11 


2.27±0.18 


M6 


17.6±2.8 


59.7±5.2 


0.59±0.09 


1.98±0.17 


Envelope** 


2334.5±32.5 




77.39±1.08 





* Measured from the spectrum 
** Aperture diameter 7'.' 62 



and the nuclear starburst is resolved into several compact 
components embedded in a larger, more diffuse envelope. 
The compact sources, in particular Ml and M2, are reminis- 
cent of young, massive, embedded star clusters. We searched 
the literature and various data archives in order to find coun- 
terparts to the MIR sources at other wavelengths, but no 
clear identifications can be made. A selection of data from 
optical to radio is shown in Fig. [5] 

The top row of Fig.Hshows two HST images, taken with 
WFPC2 and NICMOS through the F606W and F160W fil- 
ters, respectively (programmes 8597 and 7330). Both HST 
images are dominated by stellar continuum, and the ex- 
pected locations of Ml and M2 are marked with white open 
circles. No clear counterparts to Ml and M2 can be seen, 
and there is also no clear counterparts to the other fainter 
sources M3-M6. 

Although no counterparts are found in continuum op- 
tical and NIR, there is extended emission from the star- 
burst in narrow-band imag es centered on [Fe ll]1.64u.m 
and Br7 at 2.165iim JSosa-Brito et alJl20o3. fig. 31. [Fell] 
trace s shocks excite d by supernova remnants l|Vanzi fc Rickc 
Il997tl . and indeed iReunanen. Kotilainen fc Prietc1 ~ (l2003tl 
observe a high [Fen]/Br7 intensity ratio typical of shock 
excitation. This might indicate that the most massive stars 
in the starburst have already exploded as supernovae. The 
resolution in the [Fe n] and Br7 narrow-band images is how- 
ever not as good as in the VISIR image, making it difficult to 
tell where the [Fen] and Br7 emission comes from exactly. 
Some of it may come from the compact sources seen in the 
VISIR image, but it is also possible that a large contribution 
comes from the diffuse envelope. 

The M2 sour ce is clearly seen in the deconvolve d 
ll.9u.rn image by ISiebenmorgen. Kriieel fc Spoonl i2004i1 . 
but Ml is either weak, or undetected. This morphology, 
with a weak or undetected Ml and a brighter M2 source, 
is also recognized in the n arrow-band Eb (2.12u.m) image 
by ISosa-Brito et al l J200l|) . A host-galaxy subtracte d L- 
band image by fpneto^leunanen fc Kotilainen! j2002t) un- 
veils part of the circum nuclear starburst disk, probably due 
to PAH emission at 3.3 u.m. 

The bottom row in Fig. [5] shows the VISIR image and 
an ATCA radio map obtained at 8.9 GHz dMorganti et alJ 
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Figure 1. Narrow-band image and spectrum centered on the [Nell] line. The spectrum has not been continuum subtracted. The slit 
position is marked on the image. North is up and East is to the left. 
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Figure 2. Comparison with data at other wavelengths. Top row, from left to right: HST/WFPC2 F606W and HST/NICMOS F160W. 
Bottom row, from left to right: the VLT/VISIR narrow-band image and a 3cm ATCA radio map. The two circles South of the nucleus 
indicate the location of the two MIR sources, Ml and M2. The two small arrows in the HST images point to the stars that were used 
for aligning. Each major tick mark corresponds to 1 ". North is up, and East is to the left. 



dio emission may also be non-thermal synchrotron emission 
from the AGN. When the Atacama Large Millimeter Array 
(ALMA) becomes available, it will be possible to separate 
the different components at mm and sub-mm wavelengths 
and study the starburst at great detail. 



1999). Young, dust-embedded star clusters are expected to 
show thermal free-free emission at radio wavelengths, but 
the resolution in the ATCA map is ~1", not sufficient to 
match the resolution of the VISIR image. Extended emission 
from the starburst can clearly be seen, but part of the ra- 
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4 PHYSICAL PARAMETERS OF THE MIR 
SOURCES 

[Nell] is the most common Ne ion in Hn regions, ionized 
by a broad range of stellar types. It is collisionally excited, 
and the line at 12.8n.rn results from spontaneous radiative 
deexcitation. In typical star clusters, the electron density, 
n e , is a factor of 10 or more lower than the critical den- 
sity for collisional thcrmaliza tion of the 12.8n.rn trans ition, 
n c = 5.4 x 10 5 cm" 3 J Johnson fc Kobuhiickvl l2003Tl . Be- 
cause n c <C n c , the emission measure of the ionized gas can 
be estimated from the [Nell] line flux, which in turn can 
be used to predict the thermal free-free radio emission, the 
number of ionizin g photons and the number of massive stars 
jKeto et alJll999l). We proceed along the lines of the anal- 
ysis done bv lKeto et alJ dl99flf) . and calculate the emission 
measure by using that the [Nell] line intensity depends on 
the number of atoms in the upper state: 

_ g / 10 4 K\^ , hu , f n e 2 ,„„ 

1 = 3 x 10 (-zrj eM ~wJ J (^)« d ' (1) 

llKeto et alJ Il999l and references therein 2 ). The line inten- 
sity, /, has units of erg s _1 cm -2 sr~* and the emission mea- 
sure, EM — n 2 dZ, has units of cm -6 pc. We assume electron 
temperature and density of T e = 10 4 K and n e <C n c , re- 
spectively, and that the abu ndance of Ne relative to H is 
7Ne = 0.83 llKeto et al.lll999Tl . From the flux measurements 
in Table Q we derive [Nell] intensities by dividing with the 
solid angle, Q, subtended by the source, taken to be equal 
to the area of the circular aperture used for photometry. 
Table |21 lists the derived emission measures. 

The emission measure is directly related to the o ptical 
depth, t„, in the free- free continuum (iKeto et alll99sl eq. 4 
), and we estimate optical depths at 5 GHz of ~0. 01-0. 02. 
The optically thin approximation for brightness tempera- 
ture, Tb = T c t„, is therefore valid. We proceed with cal- 
culating the radio flux density at 5 GHz arising from free- 
free processes from the brightness temperature given that 
S v = 2v 2 kT h /c 2 Q,, see third column of Table H 

The Lyman continuum photon flux required to maintain 
ionization is given by 

iVphts -1 ) = 8.04 x 10 46 T- a85 [/ 3 . (2) 
The excitation parameter, U, can be written as 



U = 4.533 




flKurtz. Churchwell fc Woodlll994h . where D is the distance 
to the source. Using the radio flux calculated above, we find 
that the number of ionizing photons in the compact MIR 
sources is 2 x 10 52 s _1 as listed in the fourth column of 
Table H 

Only massive O- and B-type stars contribute to the Ly- 
man continuum flux, hence we can make an estimate of their 
numbers within each source. By making assumptions about 
the age of the compact sources and the initial mass func- 
tion (IMF), we can also derive stellar composition and mass 

2 The constant 3xl0 8 in Keto et al.'s eq. 2 is probably a typing 
error, and should be 3xl0 — 8 as used here. 



of each source. In order to do this we run a S tarburst99 
llLeitherer et aljll999l : IVazauez fc Leithererll2005T) model as 
an in stantaneous starburst with a Kroupa IMF jKroupal 
2002) with low mass cut-off at 0.10 Mq and high- mass cut- 
off at 100 Mq. Assuming age of one million years, we find 
that the compact MIR sources contain 1000-1500 O-type 
stars. The total mass in stars of each source, including the 
stars that do not contribute to the ionizing flux, is typically 
3-5 x 10 Mq. The total stellar mass (assuming an age of one 
Myr) and the number of O-stars in each source are listed in 
the last two columns of Table [5] 

The numbers vary relatively little with the age of the 
model as long as the age is < 4 Myr. As the age in- 
creases, the Lyman continuum photon flux decreases, so 
more stars are required in order to produce the observed 
photon flux of w 10 52 ' 4 s _1 . I f the average photon flux 
per OB star is 10 49 photons s~ 1 j|SmiHr J Norris fc CrowtheJ 
12001 iMartins. Schaerer fc Hillierll2005l) . the cluster cannot 
contain more than a few thousand OB stars, as derived 
above. But if the age is > 4-5 Myr, the cluster must contain 
more OB stars and hence be more massive. 

If we assume that all the ionizing flux eventually is emit- 
ted into the infrared, we can calculate an upper limit to 
the MIR con t inuum . For this, we use eqs. A2 and A3 from 
iGenzel et all <1982ft which assume that the dust emission 
can be approximated by a grey body radiating at 300 K. 
The predicted upper limits at 12.8u.m, labeled Smir in Ta- 
ble |5| are consistent with our measured 3er upper limit of 
0.27 Jy. 

The radio flux from the central regions of NGC 7582 
consists of contributions from thermal free-free emission 
from the starburst, synchrotron radiation from supernova 
remnants and the AGN. The radio flux we derive above is, 
by definition, due to thermal free-free emission and we can 
compare the predicted thermal radio emission with the ob- 
served flux to estimate the contribution from non-thermal 
processes. As existing radio measurements cannot cannot re- 
solve the different MIR components, we compare with the 
estimated thermal radio flux for the whole envelope (on the 
last line of Table^. At frequencies of 8.6, 4.9 and 1.4 GHz 
we derive free-free continuum fluxes of 33.6, 29.7 and 28.1 
mjy, respectively. The measured radio fluxes at the same 
frequencies are 49.29, 69 a nd 166 mjy llUlvestad fc Wilson! 
11984 iMoreanti et alJll999h . The contribution from free-free 
emission related to the starburst is therefore «20, 40 and 
70% in order of increasing frequency. The relative contribu- 
tion from free-free processes increases with frequency be- 
cause, as opposed to synchrotron emission, free-free pro- 
cesses produce almost flat radio spectra. We have assumed 
that there is little synchrotron contribution from super- 
novae. If the age is small enough that not many supernovae 
have exploded yet, this is a valid assumption. However, there 
may be contributions to synchrotron emission from super- 
novae in the envelope. T he observed extended [Fe 11] emission 
iSosa-Brito et al1l200lT) may indicate that some supernovae 
have already occurred. It is however possible that the com- 
pact sources are younger than the surrounding envelope. 
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Table 2. Derived quantities of the sources. The size of the envelope is taken to be equal to the area of a circle with 7'.'62 diameter. 



Source 


EM 


off 

°5GHz 


logiVph 


SWlR 


log Mtot 


logiV 




X 10 6 cm~ 6 pc 


mjy 


s- 1 


mjy 


M 


O-stars 


Ml 


2.51 


0.58 


52.39 


21.76 


5.71 


3.21 


M2 


2.44 


0.56 


52.38 


21.11 


5.69 


3.19 


M3 


1.63 


0.37 


52.20 


14.14 


5.52 


3.02 


Ml 


1.61 


0.37 


52.20 


13.93 


5.51 


3.01 


M5 


1.28 


0.29 


52.10 


11.10 


5.41 


2.92 


M6 


0.98 


0.22 


51.98 


8.49 


5.30 


2.80 


Envelope 


0.32 


29.69 


54.10 


1122.50 


7.42 


4.92 



5 CONCLUSIONS 

We have detected several compact regions of [Nen]12.8u.m 
emission in the circum nuclear starburst disk of NGC 7582. 
The compact MIR sources are embedded in a larger and 
more diffuse envelope 300-400 pc across, and do not have 
counterparts detected at continuum optical and NIR wave- 
lengths. If the sources are young star clusters, the compact 
MIR emission suggests they are embedded in their own dusty 
birth material. The MIR emission of young, embedded clus- 
ters is expected to decrease as the most massive stars ex- 
plode as supernovae and gas and dust are expelled. That 
the sources are bright in the MIR, and with no counterparts 
detected in continuum optical or NIR, suggests they may 
be young. It can however not be excluded that the clusters 
could be older and obscured by a screen of dust related to 
the circum nuclear disk. 

We estimate that the number of ionizing photons within 
each compact source is l-2.5xl0 52 s . By comparing with 
starburst models we find that the stellar mass of the two 
brighter sources Ml and M2 probably is ~5xl0 5 Mq and 
that they each contain roughly LSxlO 3 O-type stars. The 
less bright sources M3-M6 probably have smaller masses 
of 2.0-3.2xl0 5 M Q and contain ~1.0xl0 3 O-type stars. 
These properties are very similar to t hose of young, em- 
bedded star clusters fcCeto et al.lll999l : iGoriian et alJl20ol 
I Johnson fe Kobulnickv||2003ft . 

Our observation highlights the need for high-resolution 
MIR imaging and spectroscopy to study nuclear starbursts 
and individual embedded star clusters in circum- AGN envi- 
ronments. 
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